International Europhysics Conference on HEP 



PROCEEDINGS 



Discrete non-Abelian groups and asymptotically free 
models 



Sergio Caracciolo 

o 



Andrea Montanari 



E-mail: Andrea . Montanar i@lpt . ens . f r 



Andrea Pelissetto 



Recently, there has been interest in the critical behavior of two-dimensional cr-models 
in which the spins take values in some discrete subset of the sphere. In particular, two 
groups ^, ||, ||, ||] studied the nearest-neighbor er-model 

<ij> 

in which the spins have unit length and belong to the vertices of a Platonic solid, i.e. of 
a tetrahedron, cube, octahedron, icosahedron, or dodecahedron. Several quantities have 
been computed: the renormalized two-point function, the current-current correlation func- 
tion, the finite-size scaling (FSS) curve for the second-moment correlation length, and 
the four-point renormalized coupling. Surprisingly enough, the results for the icosahedral 
and the dodecahedral model are very close to the 0(3) ones, suggesting that these three 
models might have the same continuum limit. Patrascioiu and Seiler || considered 

these results as evidence for the 0(3) u-model not being asymptotically free, since the 
discrete-symmetry models have a finite (3 phase transition, which cannot be described in 
perturbation theory. However, the overwhelming evidence we have collected in the years 
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Abstract: We consider a two-dimensional cr-model with discrete icosahe- 
dral/dodecahedral symmetry. Using the perturbative renormalization group, we 
argue that this model has a different continuum limit with respect to the 0(3) o model. 
Such an argument is confirmed by a high-precision numerical simulation. 
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in favor of asymptotic freedom made Hasenfratz and Niedermayer |4], || suggest that, may 
be, the icosahedral and the dodecahedral models have an asymptotically-free continuum 
limit, in spite of the fact that the critical point is at a finite value of (3. 

Here, we wish to show that, by using some standard assumptions, the perturbative 
renormalization-group (RG) approach predicts that the suggestion of Hasenfratz and Nie- 
dermayer cannot be true. If the continuum limit of the 0(3) <r-model is correctly described 
by the perturbative RG, then any discrete-symmetry model cannot belong to the same 
universality class of the 0(3) a- model. 

The argument goes as follows |§. Consider the Hamiltonian 

H = (3 J2 ovOj-^/nK), (2) 

<i,j> i 

where <Xj is an 0(3) unit spin and I n (<Xj) is a polynomial in <Tj with the following properties: 
it has 0(3) spin n; the maxima (or minima) of I n {a x ) correspond to the set of vertices 
of a Platonic solid; it is invariant under the discrete-symmetry group of the solid. For 
all Platonic solids, it can be shown explicitly that such a polynomial exists. The model 
@l interpolates between the 0(3) (h = 0) and the discrete-symmetry model (\h\ = +oo). 
Now, with quite standard assumptions, one can show that I n { a ') is a relevant perturbation 
of the 0(3) fixed point. In other words, any arbitrarily small perturbation with discrete 
symmetry of the 0(3) cr-model drives the system to a different fixed point. 

The argument is fairly standard. Consider a p-point connected correlation function 
at zero external momenta in a finite box L 2 . If hL 2 < 1 and £ > L, we can 
compute the correlation function in perturbation theory, obtaining 

oo 

G^((3,h)=^fh^ a^(L), (3) 

i,j=0 

where t = 1/(3. The coefficients of the expansion diverge as L — > oo, since the infinite- 
volume correlation function cannot be computed directly in perturbation theory because 
of infrared divergences. However, by using the perturbative expansion (||), one can show 
that in the continuum limit G^ ((3, h) satisfies the RG equation 



a^w,h) = o, (4) 



where W(t), J^ n '(t), and y(t) are L-independent RG functions. Then, we make the follow- 
ing assumption: 

The RG equation (Q) — but not the expansion (|3|) we started from — is valid for all 
values of L, including L = oo. 

Such an assumption is routinely made in the perturbative analysis of the cr-model and is 
used, for instance, to obtain the small-i behavior of long-distance quantities, such as the 
susceptibility, correlation length, and so on. Solving Eq. (||), we obtain in the infinite- 
volume limit 

G {p) ((3,h) =G (p) (/3,0)$ (p) (z), (5) 
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where 

z = htPexp{4nrft) ~ /i£(t) 2 [log£(i)f, (6) 

p and o" are universal exponents that can be easily computed by using the perturbative 
results of Ref. @, i s a nonperturbative crossover function, and is the correlation 

length for h = 0. Precisely, Eq. (|5|) is valid in the crossover limit t — > 0, h — > 0, keeping z 
fixed. Equations (||) and (||) show that I n (cr) is a relevant perturbation with RG eigenvalue 
2, as expected on the basis of dimensional analysis. Different physical results (i.e. different 
results for universal quantities) are obtained by varying the variable z, as usual in the 
vicinity of a point perturbed by two relevant perturbations (more precisely the thermal 
direction is marginally relevant). Thus, within the standard perturbative approach, the 
discrete-symmetry model and the 0(3) model are expected to have different continuum 
limits. 

The previous argument together with the numerical results of Refs. ||, ||, || puts 
asymptotic freedom on a dangerous ground since it shows that the conventional scenario is 
wrong if the icosahedral or the dodecahedral models have the some continuum limit of the 
0(3) model. We have thus decided to extend the previous numerical work and indeed, we 
have found good evidence that the 0(3) model and the discrete-symmetry model belong to 
different universality classes: the conventional scenario is saved. However, the surprising 
fact is that these differences appear only very near to the critical point, i.e. for ^ > 10 5 ! 

In the numerical simulation we have considered the Hamiltonian (||) with 

/ 6 (<x) =u\- baj (o- 2 x + a 2 ) + 5a 2 (a 2 x + aj) 2 

+2a x a z (at - 10a 2 x a 2 y + aj) , (7) 

and h = 0.1. Such a polynomial is invariant under the rotation group of the icosahedron 
and of the dodecahedron. We measured the second-moment correlation length as defined 
m Refs. |, |, and the spin-n susceptibilities 

Xn = (*0 • *•*)) , (8) 

x 

where P n (x) is a Legendre polynomial, for n = 1, 3, 4. For each observable 0(L, j3), we con- 
sidered the so-called step function, i.e. the ratio 0(2L, ff)/0{L, /?), which, in the continuum 
limit should become a universal function of £(L, /3)/L, i.e. 

We measured the step function of the above-mentioned observables in the discrete-symmetry 
theory (i.e. keeping h = 0.1 fixed) and in the 0(3) model, thereby extending the results of 
Refs. ||, If the two models have the same continuum limit, the function computed for 
h = 0.1 and h = should coincide. 

In Figs. H|, ||we report the numerical results for the correlation length. The continuous 
line is a fit to the 0(3) data, while the points refer to the model with h = 0.1. As observed 
in previous work, there is indeed very good agreement between the numerical results for the 
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Figure 1: FSS function for the second- 
moment correlation length. 
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Figure 2: FSS function for the second- 
moment correlation length. Here, we restrict 
the horizontal range to 0.6 < x < 1.4. 
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Figure 3: FSS function for the spin-1 sus- Figure 4: FSS function for the spin-1 sus- 



ceptibility xi- 



ceptibility \i- Here, we restrict the horizon- 
tal range to 1 < x < 2. 



CD 

ho 
o 
o 



two models, but such an agreement disappears for £(L)/L > 1, where small discrepancies 
are observed. As it can be seen from Fig. |^, the icosahedral points tend to be above the 
0(3) curve and, more importantly, the discrepancy tends to increase with L: the points 
with L = 10-20 are systematically below the points with L = 50-100. 

In Figs. |l], Q the difference in behavior between the two models is quite small and not 
totally convincing. Better evidence is obtained from the results for the susceptibilities, 
since in this case the statistical errors are smaller. In Figs. |3| and [| we report the spin-1 
susceptibility and in Figs. || and ^ the spin-3 and spin-4 analogues. Again, the numerical 
results for the icosahedral and the 0(3) model agree very nicely up to £(L)/L ~ 0.8 - 
1, but then they indicate that the icosahedral FSS curve is steeper than the 0(3) one. 
Again, notice that the discrepancy between the two models increases with L, indicating 
that the observed effect is not due to corrections to scaling, i.e. it is not a lattice artifact 
disappearing in the continuum limit. 

In conclusion, the numerical results show that the icosahedral and the 0(3) model 
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Figure 5: FSS function for the spin-3 sus- Figure 6: FSS function for the spin-4 sus- 
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belong to different universality classes. Note however that discrepancies are observed only 
for £(L)/L ~ 1, which corresponds to very large values of the infinite-volume correlation 
length (see, e.g. Ref. ||). 

We thank Peter Hasenfratz and Ferenc Niedermayer for fruitful discussions. 
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